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Abstract 

The semileptonic decay form factors of the double heavy Be meson 
provide a unique opportunity to study the strong interactions between 
two heavy quarks. A fully relativistic model, with effective non-local 
quark-meson interactions, is used to compute semileptonic decay form 
factors, for both the Be and a wide range of other heavy-light mesons. 
Using these form factors predictions for decay rates and branching 
ratios are obtained. The results are compared to other theoretical 
approaches and, where available, to experimental results. In addition 
the radiative decay of B* is discussed. 
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1 Introduction 



A primary goal in the study of semileptonic decays of heavy mesons is to 
extract the values of the CKM matrix elements. The great virtue of semilep- 
tonic decays is that the effects of the strong interaction can be separated 
from the effects of the weak interaction into a set of Lorentz invariant form 
factors 1^] . Thus the theoretical problem associated with analysing semilep- 
tonic decays is essentially that of calculating the form factors. 

The focus of this work is the decay of the Be meson (for a review of the 
properties of this system see Q). This system is unique among mesons made 
up of heavy (charm or bottom) quarks, it is the only one which is stable 
with respect to strong and electromagnetic interactions. Therefore, the Be 
system is the only heavy meson for which form factors (albiet transition 
form factors rather than elastic) can be measured. These form factors then 
provide a unique probe of the dynamics of heavy quark systems. 

There are many approaches to the calculation of decay form factors, for 
example, lattice QCD QCD sum rules Q, and phenomological mod- 
elling 1^]. In this work a particular model with an effective quark- meson 
coupling is adopted. There are many models of this type ||6|, |^ |9|, p^ . 
The one used here has its genesis in the QCD version of the Nambu-Jona- 
Lasinio model |jll| extended to heavy quarks and is most closely related 
to the model used recently by Ivanov and Santorelli in a their study of 
pseudoscalar meson decays [^] . 

The advantage of this approach is that it is fully relativistic and very 
versatile. Quarks and mesons for all masses are treated within the same 
framework. For light quarks the model has the features of spontaneous 
chiral symmetry breaking and in the single heavy quark limit the form factor 
constraints of heavy quark symmetry are obtained. 

Our work differs from Ivanov and Santorelli in the choice of the quark- 
meson vertex function and in the way that parameters are fixed. A number 
of heavy mesons decays not calculated in Ref.[^] are treated here. The 
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main new results are the extension of the model to include doubly heavy 
mesons, the calculation of Be semileptonic decays and the electromagnetic 
vector to pseudoscalar transitions. 

This paper is organized as follows: the next section introduces the model, 
discusses the general method of calculation, and fits the models free param- 
eters. Sect. ^ presents the calculation of the form factors and decay rates for 
the semileptonic decays of a wide varity of heavy-light pseudoscalar mesons. 
These calculations are compared with both measured results, and other the- 
oretical approaches. Sect. ^ presents the same set of calculations for the 
eight primary semileptonic decays of the Be meson. The predictions are 
compared with other theoretical work, in order to highlight the differences 
that exist between various approaches. Sect. ^ briefly discusses the electro- 
magnetic decays V ^ P + j for a number of vector mesons, including the 
B*. Sect. ^ gives conclusions and directions for future work. 



2 Quark— Meson Coupling 

The particular quark-meson coupling used in this work is based on an effec- 
tive Lagrangian which models the interaction between mesons and quarks 
with a non-local interaction vertex |l2|. The interaction Lagrangian has 
the form 



Lint{x) = g^Mix) [dxidx2 6( X - i^ ^^^^'^^^^'^ ] 
J V \2 mi +1712 J 



f[{xi - X2)'^]qiixi)TMq2ix2), (1) 

where Tm is the Dirac matrix appropriate to the meson field M, /[(xi — X2)^] 
is a non-local vertex function, which simulates the finite size of the meson, 
and qi and q2 are the quark fields. A condition imposed on the vertex 
function is that it should render all loop diagrams UV finite. The coupling 
constant gj^.j is determined by the compositeness condition, which is the 
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requirement that the renormahzation constant of the meson fields be zero, 
i.e. 



Zm = 1 ^^^1/(^2 



= 0. 



2 dp2 

Here TlMip^) is the self energy of the meson field, given by 

[ mi-{l^ + f) 1712 -f) 

where mi and m2 are the masses of the quarks in the loop and Q is a relative 
momentum chosen to be Q = A; + ap with a = — ^ — . 

^ ^ mi+m2 

The constituent quark masses in are free parameters. As well, the 
vertex function will contain a free parameter which reflects the size of the 
meson. These parameters will be different for the different mesons. 

The use of free constituent quark propagators in expressions like @ can 
lead to a problem which reflects the lack of quark confinement in the model. 
If the meson mass Mm is greater than the sum of its constituent quark 
masses loop integrals will develop imaginary parts. This indicates a non- 
zero amplitude for the creation of a free quark-antiquark pair. There have 
been some various attempts to obviate this problem within quark-meson 



effective theories 0, @, ^, • Here we adopt the approach of Ref . |]12[ and 
use free propagators. The constituent quark masses are then fit to allow for 
the inclusion of as many mesons as possible. 

In order to carry out calculations a choice must be made for the vertex 
function /(g^). The function that was used in this analysis was the dipole 



[A2-Q2J2- 

This choice was made for two reasons; first the form of the dipole vertex 
function is the same as a propagator, allowing standard Feynman parame- 
ter techniques to be used in evaluating loop integrals. Second, the vector 
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decay constant fy would diverge if only a monopole vertex function was 
used. Since one of the primary criteria for the vertex functions is that they 
should render all diagrams UV finite, a function with UV fall-off as least as 
fast as a dipole is needed. The parameter A characterizes the finite size of 
the meson, and will be different for different mesons. To account for this 
the various values of A will be distinguished by subscripts which reflect ei- 
ther the meson type or the quark content, e.g. Ab^ and Ai^^. will be used 
interchangeably. Further, in expressions involving the vertex form factor, 
the same comvention will be used. Note that the calculations of Ref.jl^] 
used a Gaussian vertex function so that the parameters used there can not 
be compared directly with ours. 

The parameters of the model were fit to the leptonic decay constants, 
fp and fv- These quantities are defined by 

(0|-i7Vl^> = i/pP^ (4) 
{0\i^^\V,e) = M^ve'' (5) 

where My is the vector meson mass. The pseudoscalar decay constant is 
given by the one-loop expression 

[mf — {k + [771-2 ~ 

Here 7771 and 7772 refer to the masses of the quarks in the loop, this conven- 
tion will be used throughout this paper. Using the dipole vertex function, 
combining the denominators using Feynman parameters, and performing 
the integration over k yields 

f /'n-^l['^2(l - cr) +?77icr] 

fP = 9p^]Dx , (7) 

with 

a = axi + X2, 
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1] = a^xi + X2, 

A = Apxi + mlx2 + ^2X3 + (cT^ - r])Mp, 
i = k + ap, 

where Mp is the mass of the pseudoscalar meson. Likewise the expression 
for the vector decay constant is 

3At^ r ^ mims + A + a(l - ct)M^ 



= '^4^ y ^"^"^ A^ ' 

where the same defintions have been used, with the obvious change of Mp 
to My and Ap to Ay in the expression for A. 

To compute the couphng constants gp and gy, the self energies and their 
derivatives must be computed. Then the compositeness condition (^) can 
be used to find the couphngs. The self energy for a pseudoscalar meson is 
given by 

3A| f 3mim2 +p^cr(l - a) + |A 



np(p^) = ^ jDxxi ^4 — (9) 

where A = ApXi + m\x2 + rn^x^ + (cr^ — r/)p^ and all the other quantities 
are the same as the ones defined above. The self energy for a pseudoscalar 
meson is given by the tensor liy' which can be expressed as 

K^{p) = ^v{p')9,. + ^vip')^- (10) 

Unfortunately Ily 7^ fly, so this does not have the proper structure for a 
vector propagator. This problem was solved (following |jl3|) by simply drop- 
ping the riy term, which would cancel out of any calculation of a physical 
process at one-loop order (since e ■ p = 0). The relevant part of the vector 
meson self energy is given by 



^v{p') = Dxxt ^ ^ ^, (11) 
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where all the quantities appearing have been defined previously. 

The free parameters of the model are fit to the six values of fp and 
the measured values fj/^ = 0.1309 and /x = 0.075012 ]14|. These data, 



which are displayed in Table fix eight free parameters. In order to reduce 
the number of free parameters to match the available data the value of the 
strange quark mass was fixed at 450 MeV and the vertex parameter for a 
vertex containing only u and d quarks was taken (following |^]) to be 1 
GeV. In addition the following further simplifying assumptions were made 

Ams = ^ds = Ass = ^K, 
J^uc = ^dc = Asc = A/), 

Auf, = Arfb = Ksb = A-B- 

This leaves the following parameters to be fit, ruq, rric, mh, Ak, A/), A^, 
Acc, Afeb, and A^c- The parameter A^c could only be fit to a value for /b^ 
which is not in the values listed in Table |l|, hence it is retained as a free 
parameter, leaving eight to be fit. The fit to the remaining eight parameters 
is given by (all values in MeV) 

= 245, 
TUc = 1800, 
TUb = 5100, 
Ak = 1225, 
Ad = 1350, 
Ab = 1500, 
Acc = 1420, 
Abb = 2900. 

The values for the self energies, coupling constants, and leptonic decay con- 
stants arising from these parameters are displayed in Tables |^ and |^. 
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In order to fix Af,c a value of /b^ must be given. There is no experimental 
value for this quantity and theoretical estimates tend to fall in the range 



400 MeV < Jb, ^ 500 MeV (see, for example, g [T|, |T|, |Tg] ) . A further 
complication is that the mass of Mb^ is also not yet known very well. The 
current measurement |o[ is M^^^ = 6A±0.39{stat)±0.13{sys) ^ , which 
comes from the few confirmed Be events at the Tevatron. Theoretical results 
tend to lie within this range, so following the potential model prediction of 
the mass of the Be was chosen to be 6.25 GeV. 

One general argument guides the selection of Abe, it should lie between 
Ab and A;,^. With this in mind, and using the value for Mb^ above, a number 
of values of A^c were tried, spanning the possible range. Fig. |l| shows the 
value of Jbc as a function of A^c- The value selected selected for use in this 
work was A^c = 2.3 GeV, which gives /b^ = 450 MeV, a value in the middle 
of the range of the theoretical predictions. 

3 Semileptonic Decays of K, and B Mesons 

The model used in this work is phenomenological but having fixed its param- 
eters, the results for semileptonic decays are predictions. Before proceeding 
to decays of Be it is important to test the model against experimental results 
where they are available. Therefore several semileptonic decays of -ftT, D, and 
B meson are calculated. The formalism for these calculations, presented in 
this section, extends directly also to the calculation of Be decay. 

Some of the decays considered here have already been treated by Ivanov 



and Santorelli [12|. However, that work does not demonstrate the full appli- 
cablity of the approach. Apart from decays to light vector mesons, the model 
is capable of treating virtually any semileptonic decay (with the restriction 
that a value for the meson mass must be supplied as input). 
The amplitude A for a semileptonic decay is given by, 

A = ^Vqq,L,H^. (12) 
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Here Gp is the Fermi constant, Vqq' is the relevant CKM matrix element, 
is the lepton current 

and H^^ is the hadron current 

H^" = {k,e\{V^ - A^')\P), (13) 

where P is the momentum of the parent meson, k is the momentum of the 
daughter meson, and e is the polarization, if the daughter meson is a vector. 
The two currents in (|l3|) are the vector and axial A^^. If the final state 
is a pseudoscalar the hadron current can be decomposed as follows, 

{k\A^\P) = 0, 

{k\V^^\P) = U{q^){P + kr + f^{q'){P-kr, 

where f-\-{q'^) and f~{q^) are Lorentz invariant form factors. Likewise, if the 
final state is a vector meson, 

{k,e\A''\P) = f{q^)e*^' + a+{q^){e* ■P){P + kY + 

a.{q^)ie*-P)iP-kr, 
{k,e\Vf^\P) = ig{q\^^-P"el{P + k)p{P-k)^, 

where the form factors are g, /, a+, and a_. In each of these expressions 
q = (P — k) is the momentum transfer. 

For a decay to a pseudoscalar meson (with mass denoted by Mp/) the 
differential decay rate can be reduced to 

dT _ GI\Vqq>\^MI,K\^ ^ 



dq 247r3 

where, 

Mp 



IMQ'r- (14) 



K 

2 



Ml ^2 



M2 



y 



4%. (15) 
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The lepton spectrum is given by, 



da; Idtt'^ 

Cymax 

[ymax{x) - y]\f+{q^)\^) dy, 

IQ 



f 

Jo 



(16) 



with X,: 



— TTjji^- If the final state is a 

ZiVJp 



where ymax{x) 

vector meson (with mass My) the corresponding differential decay rate is 
dr _ Gl\VQQ>\''KMly 



dy 



where 



Ho 



f{q'')T2MpKg{q^), 
Mp 



(17) 



1 



M2 



y /(g2)+4i^V('z' 



and the final mass My should be subsititued for Mpi in (|15|). The expression 
for the lepton spectrum is given by 

dV _ GlI^QQ/pMl, ra(y) 



dx 



Ml 



327r3 

2(1 - 2x)[ymax{x) - y](i++{y) + 
l{y)y[2xmax -4a: + y]}, 



y+ 



(18) 



where the following definitions were made 

l{q^) = 2fiq^)g{q^), 

= ^{\fiQ')\' + MaAQ')\'-4MW\9i<l')\' 
+2{Mj.-M^-q^)f{q')a+{q^)}, 
A((?2) = {Ml-M^-q^-AM^q^, 

Xmax - 

Note that all of these expressions assume that lepton mass is zero. 
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Our model can be used to calculate all of these form factors. In all the 
following expressions the inital meson is composed of quarks with masses 
mi and m2 and the final meson is composed of quarks with masses and 
1712 (^-e. mi is the mass of the quark which decays to a new quark with mass 
ms, and m2 is the mass of the spectator). 

The form factors for decay to a pseudoscalar meson are 

^ ^ 2x 9A|,A|„ 
f+[Q ) = 9p9p' — -2 — 



3 , 



Dgx.x /^-^^'?^'--'^^ . (19) 



/ 2\ 9ApAp/ 
J-{l) = 9p9p' 2 — 



n- X- - - (72) 
Dx xiX2 ^-^^ . (20) 

The following definitions were made to simplify the expressions 
iDx = I {lldxi]6[Y,x^-?j , 

A = Ap,xi + Apa;2 + m|x3 + m^X4 + m2X5, 
_ mj 
mi + mj 

0'l,(2) = a^4,(3) +Ml2,(23)3^2,(l)) 
»7l, (2) = 3^4, (3) + /X?2,(23) ^2, (1) , 

A = A + ((7? -771 +0-1 0-2 )M|, 

+ (o-2 - ??2 + 0-iO-2)Mp, - 0-10-29^, 

K = m3{m2 - mi) + mim2 + ^{Mp + Mp, - q^) , 



e = (<7i + 0-2)0-1 Mp + (ai + o-2)o-2Mp/ — 0-10-2^^ 
Ci = mim2 



0-1 (0-1 - 1) + 0-2 ( 0-1 - 



-<72 ( 0-1 + 0-2 - ^ j Mp, + 0-2 (ai - ^ 



1 



77121713 -0-1(0-1+0-2--) Afp 



2 
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£72 (fJ2 - 1) + <Tl { £72 - - 



+0-1 ( fJ2 - - ) Q 

(€-^c)(ai±a2)±Ci+C2. 



These definitions (in addition to a and [i) will be used throughout the rest 
of this paper, with the obvious substitution of My and Ay for Mpi and Ap/ 
when the final state is a vector meson. 

The form factors for decays to vector mesons are given by 



9mp9mv 



9A^ 
7r2 



DxXlX2 



<72(m2 — ma) + ai{m2 — mi] 



777,2 



A5 



(21) 



9mp9mv {DxxiX2 X 



m\m2mz — 7(^-2 — mi — 2m3)A 



1 



(22) 



9Mp9My 



DxXlX2 



^{Pi ± p2 =f (^2m3) 
A5 



The following further definitions have been made, 
Ci = o-i(m3 - m2)(l - CTi) - miCTi, 



,1 

^3 



(72(mi - m2)(l - (T2) - m3(7f , 

m2 



1 1. 



+m3(72 I 2 ~ 

Pi = 2cri[mio-i + m2(l - fJi)], 

P2 = m2(<Ti + C72 - 1 - 2c7iC72) - miC7i(l - 2C72). 



(23) 
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Excluding the Be decays a total of sixteen pseudoscalar to pseudoscalar 
decays were considered. Due to the difficulty with confinement the cor- 
responding number of pseudoscalar to vector decays that could be treated 
was only four. Table |^ shows the predictions for the decay rates and branch- 
ing ratios for all of the decays considered. The values of the CKM matrix 
elements, and the necessary lifetimes were taken from [14|. 

Many of the decay rates treated in this section have been measured, 
hence most of the predictions can be compared to observed quantities. Ta- 
bled shows the predicted and measured results for the branching ratios. The 



experimental results are taken from [14| and the errors in the predictions 
represent the uncertainties in the CKM matrix elements. Overall, the agree- 
ment with experiment is reasonable which increases the level of confidence 
in the areas where direct comparison with experiment is not possible. 

Table ^ shows values of /+(0) as computed in this work and in various 
other theoretical approaches. The other approaches are widely varied: [^] 
uses the ISGW model, uses the WBS model, |Q gives results from a 
bag model, [^] uses a Dyson-Schwinger equation approach and [Q] gives 
lattice QCD results. Of particular interest is the work of Ivanov et al. [^] 
which uses the quark confinement model. This quark-meson model is based 
on similar considerations to the model used in this work so its predictions 
should be close to ours. 

As well, the decay B D+i^+i^i can be treated in a model independent 



way using the HQET. Ivanov et al. have shown in several papers [ffO, 12 



that quark-meson models of the type used here give the correct tree level 
HQET relations in the infinite mass limit. Nevertheless a direct check with 



finite quark mass is useful. The HQET gives the prediction |26| 
fH^ET 2^ = = L138, 

which compares well with our value f+{qfnax) — 1-133. 

The most important comparison that can be made is with Ref.|12]. This 
paper uses a different vertex function to treat B and D decays. This serves 
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as a check on the dependence of the model on the choice of vertex func- 



tion. Apart from the case B tt + £ + u agreement with |12] is very 



good. In addition [12| presents the values of f+{q'^) over the full range of 
q'^. Overall agreement is good between the two calculations, Fig. |2| illus- 
trates the agreement in the case D° ^ K- + e+ + Uf,. Fig. I shows the case 
TT~ -\- i"^ + U£, for which the agreement is better over the whole range 
than indicated in Table |5[ 

Due to lack of confinement very few pseudoscalar to vector decays can 
be calculated. Of the few decays treated in this work only the decay B — > 
D* + £ + i>£ has been studied extensively. Table ^ compares our predictions 
with some other calculations. Overall the agreement is reasonable. 

4 Semileptonic Decays of the Be Meson 

The methods of the previous section can be directly applied to the semilep- 
tonic Be decays. Using the procdeure outline above, decay rates, lepton 
spectra, and branching ratios can be computed. In this work the lifetime 
of the Be was taken to be 0.5 ps, which agrees with the CDF value of 
T^DF ^ OAGtoll ± 0.03 ps |0[. Table shows /+(0), the total 



decay rate F and the branching ratio for the four pseudoscalar decays. For 
the decays to vector mesons, values of the form factors at q'^ = as well as 
total decay rates and branching ratios are displayed in Tables ^ and |9[ 

There are a number of other calculations of the semileptonic decays of 
Be- A comparison of some results for the dominant decay modes is given 



in Table 10. In contrast to the situation in Sect. |^ where our quark- meson 
model predictions, for the most part, agreed with other models and the 
various other models agreed with each other, there are substantial differences 
between calculations of Be decays. The clearest examples of this are the 
predictions for the decays to the B* and J/'0. These two decays are expected 
to be the most important semileptonic decay channels However there is 
disagreement not only over the values of the branching ratios but also as to 
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which decay wih be favoured. For example, the quark-meson model used 
in this work predicts the decay to the J/ij) to be slightly favoured over the 
decay to the i?* while the heavy quark approach used in and [^] 
predicts the decay to i?* to dominate. This divergence of predictions may 
be expected; the heavy-heavy quark content of the Be poses a challenge 
for models. Light-quark mesons may be constrained by chiral symmetry 
and heavy-light mesons by heavy quark symmetry. On the other hand 
the physics of heavy-heavy systems is less constrained by symmetries so 
extending models into this domain provides a severe test. 

5 Electromagnetic Decays V ^ P + ^ 

In addition to semileptonic decays the electromagnetic decays of vector 
mesons can be treated within our effective quark-meson coupling model. 
Since the amplitude involves the matrix element it is clear this 

process will be related to the form factor g{q^). The the amplitude for this 
process is 

A = -2ee^''"'l'~e^e,py^ppp[Qigi{Q) + Q252(0)], (24) 

where Qi (2) is the charge of qi^ (2); and pp y are the momenta of the pseu- 
doscalar and vector mesons. The functions 5i(0) are the form factors given 
by (^11), with the appropriate masses inserted, and with = 0. The appro- 
priate masses in these functions are given by the interchange of Mp and My 
and the subscript which denotes which of the quark lines the gauge field is 
coupled to {i.e. for 51(0) the appropriate expression sets m3=mi). Defin- 
ing Qyp^ = 2[Qigi + (5252] 1 and summing over initial and final polarizations 
gives 

2 

9VPT (25) 



2a'K 



Ml 
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where a = -Xf is the fine structure constant. Standard techniques 1 31 1 yield 



137 

the total rate ^ 

^yp^ = YA^y9lp,[^-^) ■ (26) 

Electromagetic decays have been the subject of several theoretical stud- 
ies. As well the decay J /if) rjc + ^ has been measured. Table ^ shows 
our predictions for gyp^ along with the single experimental result and the 
predictions of some other models. In |32] and two different heavy quark 



approaches were used. The quark confinement model which has some 
similarity to the quark-meson model used in this work, gives results which 
are quite close to ours. 

There are measured branching ratios for the D* decays, however no 
lifetime measurement has been made. Therefore our predictions (which do 
not include the lifetime) cannot be compared directly with experiment. In 
order to obtain branching ratios a theoretical estimate of the lifetime must be 
used. The quark confinement model is the ideal choice, since its predictions 



are closest to our work. Using the results from [13| and our predictions 
for the total rates (obtained from (|2^)) the following branching ratios are 
obtained: 

BR[{D*f ^ D° + j] = 33.0%, 
BR [{D*)+ + j] = 1.43%. 

These compare well with the experimental values IT 



0, 



BRe,cpt.[{D*f ^ D^ + j] = 38.1% 
B Re.pt. [{DT ^D+ + j] = 1.1%. 

In order to treat the electromagnetic decay B* Be + 'j the mass of 
the B* meson must be specified. Theoretical estimates indicate that the 
mass difference should be small; Mb* — Mb^ < 100 MeV. To examine the 
effect of a small change in the B* mass, the self energy and coupling constant 
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were calculated over a range of masses. These results were used to calculate 
dB'Bc'y ^'^'-^ displayed in Table 12. The decay rate is shown in Fig. ^ 



The radiative decay of B* has not been studied extensively. A QCD sum 



rule approach |34], using Mb* = 6.6 GeV and Mb^ = 6.3 GeV, gives the 
result Qb^Bc"! ~ 0-270 ± 0.095Gey~^. Using these masses our prediction is 
dB'Bc-y ~ 0.2196Gey~^. The two values are in agreement. 

6 Conclusion 

A Lagrangian which models mesons in terms of an effective non-local quark- 



meson interaction vertex 12 1 was extended in this paper to describe 
mesons such as B* composed of two heavy quarks. The model has the 
advantage of treating all quarks (heavy and light) on the same footing, 
thereby permitting a unified investigation of heavy — > heavy, and heavy 
— > light quark decays. The model does not provide a complete dynamical 
description of quark interactions, meson masses can not be calculated and 
must be introduced as input parameters. Quark masses and the parame- 
ters associated with the quark-meson vertex were determined by fitting the 
pseudoscalar and vector meson decay constants fp and fv- Due to lack of 
confinement in this approach some light vector mesons had to be excluded 
from the analysis. 

To both test the model, and demonstrate its versatility, a large number 
of semileptonic decays of K, D, and B mesons were analysed. Agreement 
with measured results and other theoretical approaches was good. 

The main focus of this work was on the analysis of the semileptonic 
decays of the doubly-heavy Be meson. The form factors characterizing the 
strong interactions of the Be system were computed over the entire available 
range of momentum transfer. Using these results, decay rates and branching 
ratios were computed for all eight decay going to mesonic ground states. 

A comparison with some other approaches highlighted the significant 
differences among various model predictions concerning the Be- On very 
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general grounds one would expect the decays to J /if) and {B*)^ states to 
be the most important. However, there is no agreement from models which 
of the channels dominates and absolute rates differ by a factor of 3 to 4. 
This is in contrast to the situation in K, D and B meson decays where much 
smaller differences between different models are found. 

As a further illustration of the versatility of the model, the electromag- 
netic decays V ^ P + ^ were investigated. A reasonable description of 
J/ "0 and D* radiative decay was found and the rate for B* ^ Be + j was 
calculated for a range of B* masses. 

There is room for further analysis within the model considered in this 
work. Hadronic decays, such as B^ J /ip + vr'^ can also be treated. 
A detailed analysis of all of these decays, combined with the results for 
the leptonic and semileptonic decays, could be used to make a prediction 
for the lifetime r^^. A further area that needs work is the difficulty with 
confinement. The natural solution to this problem appears to be provided 
by the quark confinement model Recently Ivanov et al. have proposed 
a modification of the quark confinement model which may aid its application 



to heavy mesons |35|. 
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Table 1: Pseudoscalar Decay Constants 



Meson (P) 


fp (MeV) 


Source 


7r+ 
K+ 
D+ 

Dt 
B+ 


130.7 ±0.5 

159.8 ± 1.8 

loq+12+41+9 
-'■O0_^3_Q_25 

090+10+51+3 
^^y_ll_0-19 

1 c;r+12+29+9 
-•■00-14-0-9 
-, 77+11+39+13 
'-12-0-11 


Expt. [14| 
Expt. [14| 
Lattice QCD H 
Lattice QCD 1T| 
Lattice QCD |l§ 
Lattice QCD H 



Table 2: Fitted Values of Pseudoscalar Meson Properties 



Meson 


n(M2) (GeV2) 


9p 


fp (MeV) 


Tr+ 


0.038588 


5.14916 


131.06 


K+ 


0.068760 


5.16139 


160.85 


D+ 


0.083122 


6.26671 


182.80 


Dt 


0.087025 


6.82039 


223.81 


B+ 


0.134842 


6.44510 


142.21 




0.136072 


7.75408 


187.42 



Table 3: Fitted Values of Vector Meson Properties 



Meson 


n(M2) (GeV2) 


9v 


fv 


T 


0.095419 
0.621464 


8.81505 
7.54393 


0.131227 
0.074796 
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Table 4: Predictions for Decay Rates and Branching Ratios 



Parent 


Daughter 


r (ps 


) 




BR{%) 


on-exptX/o) 


K+ 




3.589 X 


10' 


f: 

-u 


4.45 




3.18 ± 0.08 




vr" 


7.274 X 


10- 


-D 


0.0650 




0.0670 ± 0.0007 




7r~ 


7.274 X 


10- 


-D 


37.6 




38.78 ± 0.27 




7r~ 


5.488 X 


10- 


Q 


0.228 




0.37 ± 0.06 




K- 


8.476 X 


10- 


-z 


3.52 




3.50 ± 0.17 


D+ 


n 


2.790 X 


10- 


q 


0.295 




0.31 ± 0.15 


D+ 




8.515 X 


10- 


-2 


9.000 




6.8 ±0.8 






4.184 X 


10- 


_3 


0.195 






Dt 


Z)0 


4.786 X 


10" 


-8 


2.24 X 10" 


-6 






vr" 


6.455 X 


10^ 


~5 


0.0101 




0.018 ±0.006 




D- 


1.716 X 


IQ- 


-2 


2.68 




2.00 ± 0.006 


B° 


(D*)- 


3.983 X 


10" 


-2 


6.21 




4.60 ±0.27 


B+ 


tt" 


3.457 X 


10" 


-5 


5.70 X 10" 


-3 


< 0.22 


B+ 


DO 


1.726 X 


10- 


-2 


2.85 




1.86 ±0.33 


B+ 


(Z)*)0 


4.059 X 


10" 


-2 


6.70 




5.3 ±0.8 


B° 


K- 


6.118 X 


10- 


-5 


9.42 X 10" 


-3 






D7 


1.642 X 


10- 


-2 


2.53 






B's 




4.185 X 


10- 


-2 


6.45 






B's 


B- 


2.619 X 


10- 


-8 


4.03 X 10" 


-6 




B's 




8.032 X 10- 


10 


1.24 X 10" 


-7 
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Table 5: Comparison of this work with other approaches for the form factor 
/+(0). Here P is the parent meson and D is the Daughter 



p 


D 


This work 


K+ 




0.9617 




K~ 


0.7869 






0.6292 




D- 


0.7977 






0.2848 


B° 


K~ 


0.2452 







i 



0.76 
0.69 
0.69 
0.33 



0.780 
0.681 
0.684 
0.293 



0.71 
0.8 

0.33 
0.36 



0.78 



0.53 



0.98 
0.74 

0.73 
0.51 



0.73 
0.65 

0.27 



Table 6: Vector Form Factors for the Decay B —>■ D* + + f £ 



Reference 


5(0) (GeV-^) 


a+(0) (GeV-^) 


/(O) {GeV) 


This Work 

i 

|27| 
[23| 


-0.10391 
-0.09745 

-0.16 
-0.09054 


-0.09240 
-0.09471 

-0.15 
-0.07271 


5.286 
4.736 
6.863 
3.863 



Table 7: Predictions for Be ^ P Decays 



p 






r (ps-^) 


BR (%) 


B^ 


0.4504 


0.6816 


9.7001 X 10- 


-4 


0.049 


B^ 


0.5917 


0.8075 


1.8774 X 10" 


-2 


0.94 




0.1446 


1.017 


2.8244 X 10" 


-5 


0.0014 


Vc 


0.5359 


1.034 


1.0355 X 10" 


-2 


0.52 
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Table 8: Predictions for Be ^ V Decays 



V 


r (ps-^) 


BR(%) 




1.048 X 10~3 


0.052 




2.872 X 10^2 


1.44 


{D*f 


4.739 X 10^^ 


0.0024 




2.943 X 10^2 


1.47 



Table 9: Predictions for the form factors at = for B^ ^ V decays. 



V 


5(0) GeV-^ 


a+(0) GeV-^ 


/(O) GeV 


{B*f 


-0.1671 


-0.0463 


3.383 


(Btr 


-0.2402 


-0.0673 


5.506 


(L>*)0 


-0.0211 


-0.0127 


0.8296 


J/i; 


-0.0784 


-0.0543 


4.918 



Table 10: Branching Ratios for the Semileptonic Decays of the Be- 



Decay Meson 


This Work 


1 


i 




28| 


1 


m 




3C| 


B's 

(Btr 

Vc 


0.94% 
1.44% 
0.52% 
1.47% 


(0.68 ± 0.23)% 
(0.68 ± 0.23)% 
(0.57 ±0.17)% 
(0.68 ±0.17)% 


1.35% 
3.22% 
0.553% 
2.32% 


0.836% 
2.13% 


0.80% 
2.3% 

0.15% 
1.5% 
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Table 11: The decay constant gyp-y in GeV 



V 


This Work 


Expt. 


14 


1 


1 


23 




1 


32 


1 




33] 




13 




{D*f 
{D*)+ 

Jl^ 
{B*f 
{B*)+ 

(Btr 

T 


2.0321 
0.5224 
0.2369 
0.7419 
0.9770 
1.7627 
0.6417 
0.1314 


0.5538 


1.043 
0.1535 

0.3098 
0.3461 


1.079 
0.0702 

0.4177 
0.6540 


1.598 
0.2886 
0.1917 

0.5720 
0.9701 


1.939 
0.3950 
0.2598 

0.9104 
1.618 



Table 12: Calculations of B* Properties and gB*Bc 



Mb* (GeV) 


n((MBj)2) {GeV') 


OB* 


9B*B.y (GeV-') 


6.25 


0.2843 


8.627 


0.3013 


6.26 


0.2877 


8.541 


0.3011 


6.27 


0.2912 


8.455 


0.3008 


6.28 


0.2947 


8.369 


0.3006 


6.29 


0.2984 


8.283 


0.3003 


6.30 


0.3021 


8.196 


0.3001 


6.31 


0.3059 


8.109 


0.2999 


6.32 


0.3097 


8.022 


0.2995 


6.33 


0.3137 


7.934 


0.2993 


6.34 


0.3178 


7.846 


0.2990 


6.35 


0.3220 


7.758 


0.2987 
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Figure 1: /b^ vs. A^c 
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Figure 4: The decay rate T [B* B^ + j] vs. M^.. 
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